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Abstract. In this paper, we give a new type of a posteriori error estimators suitable 
for moving finite element methods under anisotropic meshes for general second-order 
elliptic problems. The computation of estimators is simple once corresponding Hessian 
matrix is recovered. Wonderful efficiency indices are shown in numerical experiments. 
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1 Introduction 

Nowadays adaptive algorithms have been an indispensable tool for most finite element 
simulations. They basically consist of the ingredients "Solve - Estimate error - Refine 
mesh" which are repeated until the desired accuracy is achieved. Generally, they can be 
classified into three types: h-, r- and /ip-version. In this paper we consider the second in- 
gredient (Estimate error) for r-version(or moving finite element method) under anisotropic 
meshes. 

Then, what does "anisotropic mesh" mean? Denote by the diameter of the finite 
element K, and by qk the supremum of the diameters of all balls contained in K. It is 
assumed in the classical finite element theory that 



(The notation < means smaller than up to a constant.) Elements which satisfy (jl.ip are 
called isotropic elements. 

Many physical problems exhibit a common anisotropic feature that their solutions 
change more significantly in one direction than the others. Examples include those hav- 
ing boundary layers, shock waves, interfaces, and edge singularities, etc.. In such cases 

*Department of Mathematics, China University of Mining and Technology, Beijing 100083, China 
(email: yinxbQlsec. cc.ac.cn) 

■^LSEC, ICMSEC, Academy of Mathematics and Systems Science, CAS, Beijing 100080, China (email: 
hhxieQlsec. cc.ac.cn) 



hK < QK- 




it is advantageous to reflect this anisotropy in the discretization by using meshes with 
anisotropic elements (sometimes also caUed elongated elements). These elements have a 
small mesh size in the direction of the rapid variation of the solution and a larger mesh 
size in the perpendicular direction. That is to say, anisotropic elements do not satisfy 
condition (jl.ip . Conversely they are characterized by 



where the limit can be considered as /i — t- (near edges) or e — > (in layers) where 
e is some (small perturbation) parameter of the problem. Indeed anisotropic meshes 
have been used successfully in many areas, for example in singular perturbation and flow 
problems [21 El E HI [Ml EI] and in adaptive procedures [9l[l2l[2ll[2Z]- For problems with 
very different length scales in different spatial directions, long and thin triangles turn out 
to be better choices than shape regular ones if they are properly used. This motivated an 
intensive study on the error analysis for anisotropic meshes in the finite element method. 
For instance, Apel [3] described an error estimate in terms of the length scales hi and ^2 
along the x and y direction, respectively. Berzins [8] developed a mesh quality indicator 
measuring the correlation between the anisotropic features of the mesh and those of the 
solutions. Kunert |2U) introduced the concept of "matching function" which measures the 
correspondence between an anisotropic mesh and a given function. Using this concept he 
gave three types of a posteriori error estimator for anisotropic meshes under the assumption 
that the anisotropic mesh is 'adapted' to the anisotropic solution. Formaggia and 
Perotto [13] used the spectral properties of the affine mapping from a reference triangle 
to obtain a full information about the orientation, dimension and aspect ratio of a given 
element. After that they proposed a posteriori estimators for elliptic problems under 
anisotropic meshes. Picasso |25j combined the method in |15j and a Zienkiewicz-Zhu error 
estimator to approach the error gradient. Cao [11] revealed the precice relation between 
the error of linear interpolation on a general triangle and the geometric characters of the 
triangle. This list is certainly incomplete, but from the papers we can find the interpolation 
error depends on the solution and the size and shape of the elements in the mesh. 

In the mesh generation community, the error estimate is often studied for the model 
problem of interpolating quadratic functions. This model is a reasonable simplification of 
the cases involving general functions, since quadratic functions are the leading terms in 
the local expansion of the linear interpolation errors. For instance, Nadler [23] derived 
an exact expression for the L^-norm of the linear interpolation error in terms of the three 
sides £1, £2, and £3 of the triangle K, 



where |-ftr| is the area of the triangle, di = £{ ■ H£i with H being the Hessian matrix of u. 
Bank and Smith [5| gave a formula for the if^-seminorm of the linear interpolation error 
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where d = [di, ^2,(^3]^, 

^=48L^ 24 -£2 \h? + \^2? + \h? 2^2-^3 

' ' V 2^2-4 |4p + K2p + |4|' / 

In this paper we'll develop the formula for iJ^-seminorm of the linear interpolation error 
l|V(n-.x,)||i.(^)« ^^Y^cln\ (1.5) 

where Cj = ^j+i • H£i^2, and that of discretization error 

1 ^ 

l|V(u-«/,)||i.(^) « E E {fK+mdnUhk)di. (1.6) 

The quality of an a posteriori error estimator is often measured by its efficiency index, 
i.e., the ratio of the true error and the estimated error (in some norm). An error estimator 
is called efficient if its efficiency index together with its inverse remain bounded for all 
mesh-sizes. It is called asymptotically exact if its efficiency index tends to one when 
the mesh-size converges to zero. Prom numerical results we see our estimators are often 
asymptotically exact although we couldn't prove it rigorously. 

The paper is organized as follows. In section 2 we give some preliminary results, espe- 
cially the error expansions for u — uj and \/{u — uj). In section 3 these error expansions are 
used to derive a posteriori error estimators for the interpolation error and discretization 
error, respectively. Section 4 contains "efficient index" tables and pictures from numerical 
experiments for some second-order elliptic problems which yield anisotropic solutions. The 
results show remarkable agreement with the theoretical predictions. Finally, in section 5 
we state our conclusions and direction for further research. 



2 Preliminaries 

Consider the following model problem. Find u: (Z7i? ^TZ such that 

, Lu = -E^£:{av^)+bu = f inO, 
u = on dfl, 

where b = 6(x) > a.e. in Q and aij = aij(x) are given functions. The domain Q 
is an open, bounded subset of TZ'^ and the operator L is elliptic and self-adjoint. The 
corresponding variational formulation seeks u E Hq (Q) such that 

a{u,v) = {f,v) \fveV = H^{n), (2.2) 
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where 



a 




+ buy I dx 



and 



if,v)^ [ fvd: 



We shall use the standard notations in [13] for the Sobolev spaces H^{Q) and their asso- 
ciated inner products {•,-)s, norms || • \ \s, and seminorms \ ■ \s for s > 0. 

By = {Th} we denote a family of triangulations Th oi Q,. Let be the space of 
continuous, piecewise linear functions over Th, and Vb,h = Vh ni?o(r2). The finite element 
approximation problem of ()2.2p seeks Uh € Vo,h such that 



The three vertices of an arbitrary triangle K ^ Th are denoted by ai = a2 = 

(3^2,2/2)"^ and as = {x3,ys)^. Additionally we define the edge vectors ii = sls — a2, 
^2 = ai — as and ^3 = a2 — ai (Figure 1). Denote by uj the linear interpolation of u at 



the three vertices of K. Let {Aj(x)}?^^ be the barycentric coordinates of K. From [131I26J 

we know for a quadratic function u over K the following formulas hold: 



a{uh,Vh) = if,Vh) ^Vh^Vo^h- 



(2.3) 




Figure 1: notations in a single element K. 




u(x) — it/(x) 



- ^ Ai(x) [(x - a;) • H{x - ai)J Vx G K 



(2.4) 



1=1 



V( n(x) — u/(x) 



) 




(2.5) 



i=l 



where H is the Hessian matrix of u. 
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3 A posteriori error estimates 



Many authors have discussed the interpolation error to derive their adaptive algorithm ( 
[H[5|[8|[T5tll9j). However, the interpolation error is different from the discretization error 
in most cases. In this section we first discuss the former and then the latter. Finally we 
will analyze their relationship using the concept "superapproximation" . 



3.1 An a posteriori error estimator for the interpolation error 



Theorem 3.1. Letu be a quadratic function anduj is the Lagrangian linear finite element 
interpolation ofu. Denote by H the Hessian matrix of u. The following relationship holds: 

1 ^ 

l|V(n-n,)||i.(^) = — — J^cfK, 



|2 



(3.1) 



i=l 



where Ci = ii+i ■ H£i+2- Here we prescribe i + 3 = i,i — 3 = i. 
Proof. From (j2.5p . we have 

VAi(x) (x - Ri) ■ H{^ - aj 



\V{u-Ul)\\^2(^K) 



K 



dx. 



1=1 



Due to the properties of the barycentric coordinates it is known that 
3 



^V\(x) 



fx - ai) ■ H(x - a,; 



i=l 



G P2{K). 



(3.2) 



(3.3) 



We use the second-order quadrature scheme which is exact for polynomial of degree less 
or equal to 2, i.e. 



K 



(f{x)dx = ^ K (f{a.i2) + (p{a.23) + ^^(asi) , Vv? G P2{K), 



(3.4) 



Vi+l - yi+2 Xi+i - Xi+2 



2\K\ 



2\K\ 



where ajj is the midpoint of the segment a^aj. Notice that 
VAi(x) = 

after a simple calculation we get 
Here we set 

m = 



1 



48|K| 
KeTh I I i=i 



□ 



(3.5) 



as the a posteriori estimator for ||V(n — ui)\\i^2(yiy 

Remark 1 Using the same technique we can also get the corresponding estimator for 
\\u-ui\\q^q, denoted by ry/o- 

Remark 2 The estimators r// and r/70 have been given in different forms, for example, 
in ElflllEa], e.t.c. 
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3.2 An a posteriori error estimator for the discretization error 



In this subsection an a posteriori error estimator for the discretization error of problem 
([^ will be given. 

Theorem 3.2. Assume L be an elliptic and adjoint operator, c(x) and aij{x) be zero and 
constant functions, respectively. We have the following estimate: 



3 

•^^ Ki^Th i=l 

where fx = //^ /(x)dx and [dnUh]ii is the jump of the conormal derivative of Uh 



(3.6) 



dnUh = 



duh 



across the edge £{, with n = (ni,n2)"^ the unit outward normal vector. 
Proof. Using the Galerkin orthogonality, we have 



V(n — Uh)\ (ix = a(n — Uh, u — Uh) = a{u — Uh,u — uj) 



a,; 



d{u-Uh) d 



Ken 



{u - ui) 



dx 



d ( d{u-Uh) 



K L — dxi V''^ dxj 



(u - Ul) 



dx+ I dn{u - Uh){u - uj)ds\ 

dK ^ 



fx 



^{/ [ i^nUh]{u - ui)ds^ + "Y I (^f - jj^j{u-ui)dyL 

I / T^(^-"/)^x + ^ / [dnUh]{u - ui)ds\ 
77^^ Wir 1-^1 ^ JdK J 

^33 ^3 

U E {fKY.d^ + Y.\^^\[^nUh],,d^ = -- ^ J] (/i^ + | M,,) 



KeTh 
1 

~24 



K&Th i=l 1=1 



where we use the error expansion (j2.4|) and the second-order quadrature scheme on K and 
dK, respectively. □ 

3.3 Discussion of the estimators 

From the Theorem 13.21 we get easily an a posteriori error estimator for the discretization 
error: 



1 



K&h i=l 
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Obviously this estimator can be computed easily provided that H is properly given. A 
number of numerical recovery approaches have been proposed in the literature for second- 
order derivatives p!t l221l29p30p 32j . Comparisons of these techniques have also been made in 
[T0|l28] . Particularly the authors [28] compared four methods for reconstructing the second- 
order derivatives of a piecewise linear function: DLF(DoubIe linear fitting), SLF(Simple 
linear fitting), QF(Quadratic fitting) and DL2P(Double L^-projection). In this paper we 
will recover H using the quadratic fitting method elaborated by Zhang [29]. 

To end this section, it is advantageous to discuss the relationship between the interpo- 
lation error \\V{u — ui)\\o^n and the discretization error ||V(u — m/i)||o,q- 

Denote by N the number of elements in Th- Assume ||V(n — u/)||o,q ~ CN'^/'^ and 
||V(ti — ti/i)||o,f2 ~ CN'^/'^. Then, by simple calculus, we have 

\\^{u-ui)\\l^-\\V{u-Uh)\\l^= [ V{ui + Uh-2u)-{Vui-Vuh)dx. (3.7) 



From (|3.7p we conclude that if 

\\V{uj-Uh)\\o,n<CN-i-\ (3.8) 
where 7 > 0(this phenomena is called super approximation [SlEI]), then 

l|V(n - ui)\\l^ - ||V(n - UH)\\ln = 0{N-'-^). (3.9) 
Assume rjj be an asymptotically exact estimator of |lV(ti — ti/)||o,Q, that is 

lim T— — — — — — — = 1. 

iv^-oo ||V(n - ui)\\o,n 

Then 7]j can also be used as the estimator of || V(u — n/j)||o,n because 

lim - — ; — -7^ — = lim - — —7^ ; = 1. (3.10) 

N-^oo\\V{u-Uh)\\l^ l|V(u-n/)||2^f^ + 0(iV-i-^) ^ 

However, the superapproximation can be proved only in some structured meshes such as 
uniform and uniform Chevron triangular meshes in [21], and 0{h'^") irregular triangular 
meshes in [6], under the assumption that u is very smooth. When the solution doesn't 
have superapproximation property we couldn't replace the discretization error by the 
interpolation error. Fortunately, from numerical experiments in section 4 we guess that 
the superapproximation always holds during the adaptive procedure. 

3.4 Problem for general coefRcients 

For the discussion above we assume that 6(x) and aij(x) are zero and constant functions, 
respectively. In fact, we can get the corresponding results for the general smooth functions 
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6(x) and ajj(x), if we use the higher order quadrature scheme and notice that 



d{u - Uh) d 
dxj dxi 



(u - ui) 



dx 



^^Tfi i,j — 1 



dx+ I dn{u — Uh){u — uj)ds> 

dK ^ 



+ 



Ken 



_d_ 



+ 



dK 



{dn^ - dn){u - Uh){u - ui)ds^ 



d{u - Uh] 



dxj 



{u - Uj) 



dx 



dx+ I dn{u — Uh){u — uj)ds> , 
dK J 



where aL{x) is the Lagrangian hnear finite element interpolant of ajj(x), and 



1 d 
x) = -r^ / aij{x)dx, dnUh = ^ aij{x)—-^ni, dn^Uh = ^ aij(x) 

' ' "^-^ i,i=l i,i=l 



dx,, 



-Ui 



4 Numerical experiments 



First we give some definitions, 
E=.. . , EI- 



\^{'^-Uh)\\ln 



. ^ — , Where H is exact Hessian matrix, 
\V{u-UhY^ 



-, EIr 



2 

m 



-, Where Hr is recovered by Zhang 



l|V(n-n,^)||2_^' \\V{u-Uh)\\l^'' 

We induce the exact Hessian for comparison in examples 4.1-4.4 and 4.6, while in example 
4.5 we just use the recovered Hr where the true solution doesn't belong to //^(Jl). 

Because we use the Hessian recovery technique in [29J, it is advantageous to show how 
this technique works. From this point we will verify if there exists a positive number d 
such that \\H — Hr\\o,n < CN^^ in our numerical experiment. 

Example 4.1 This example is to solve the boundary value problem of Poisson's equation 



Au 



/, xGl^=(0,l)x(0,l), 



(4.1) 



where the Dirichlet boundary condition and the right-hand side term are chosen such that 
the exact solution is given by 



U{x) = ( 1 -I- e 2e j 



(4.2) 
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with e being taken to be 0.005(taken from [19j). Here we use the Delauney mesh generator 
to get the nearly uniform mesh, where nu is the number of initial points on the boundary. 
See Table 1 and Figure 2 for more details. 




1 1.5 2 2.5 3 3.5 4 4.5 5 




(b) 



Figure 2: (a) Four estimators and (b) the initial mesh of the example 4.1. 



Table 1: Four estimators and S in example 4.1 



nu 


N 


E 




EI 


EIr 


\\H -Hr\\ 


5 


40 


3744 


0.395664 


0.361983 


0.475616 


0.250096 


172.773 




80 


8664 


0.701157 


0.584010 


0.621573 


0.434262 


92.8695 


1.48 


160 


15154 


0.799865 


0.701711 


0.726035 


0.560006 


62.6707 


1.41 


320 


23674 


0.875311 


0.795383 


0.815560 


0.680712 


40.1777 


1.99 


640 


34108 


0.910547 


0.845726 


0.862312 


0.746308 


29.0403 


1.78 



From example 4.2 to 4.6 meshes are generated using a C++ code BAMG(Bidimensional 
Anisotropic Mesh Generator) developed by Hecht [17] (r— version adaptive procedure). 
Example 4.2 The same problem as in example 4.1. In fact the solution exhibits a sharp 
layer on line xi + X2 — 0.85 = 0. The result is shown in the following table, where p stands 
for the step of the adaptive procedure. Results are listed in Table 2 and Figure 3. 



Table 2: Four estimators and 5 in example 4.2 



p 


N 


E 




EI 


EIr 


\\H-Hr\\ 


S 


1 


94 


-0.015593 


0.136566 


0.812635 


0.127626 


265.946 




2 


113 


0.035375 


0.176612 


0.849958 


0.179425 


233.315 


1.42 


3 


189 


0.668308 


0.540584 


0.601845 


0.414778 


153.646 


1.62 


4 


272 


0.950449 


0.846033 


0.929412 


0.738008 


62.0580 


4.98 


5 


278 


0.994030 


0.916821 


1.018023 


0.850781 


34.0592 


55.0 



Example 4.3 This example is to solve the boundary value problem of Poisson's equation 

-Au = f, x G S7 = (0,1) X (0, 1), (4.3) 
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1 1.5 2 2.5 3 3.5 4 4.5 5 

(a) (b) 




Figure 3: (a) Four estimators, (b) the final mesh and (c) Uh of the example 4.2. 

where the Dirichlet boundary condition and the right-hand side term are chosen such that 
the exact solution is given by 

u(x) = e^?-°-^ (4.4) 

This is an extreme example for anisotropic behavior where the function u only depends 
on one variable xi or X2- Such functions are the real challenge in the a posteriori error 
analysis since one is not allowed to use this knowledge. It is obvious our four estimators 
perform very well. See Table 3 and Figure 4 for more details. 



Table 3: Four estimators and S in example 4.3 



step 


N 


E 




EI 


EIr 


\\H-Hr\\ 


5 


1 


26 


0.863888 


0.659135 


1.14815 


0.700543 


1.02272 




2 


26 


0.899800 


0.695822 


1.13359 


0.721021 


0.976705 




3 


32 


0.984515 


0.804626 


1.13656 


0.835906 


0.808743 


1.82 


4 


43 


0.993768 


0.938219 


1.08240 


0.967953 


0.402943 


4.72 


5 


66 


0.990844 


0.960663 


1.04390 


0.970201 


0.190728 


3.49 
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(a) 



(b) 
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0.9 
0.8 
0.7 
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0.5 
0.4 



Figure 4: (a) Four estimators, (b) the final mesh and (c) Uh of the example 4.3. 

Example 4.4 This example is to solve the boundary value problem of Poisson's equation 

-Au = f, XG17 = (-1,1) X (-1,1), (4.5) 

where the Dirichlet boundary condition and the right-hand side term are chosen such that 
the exact solution is given by 

u(x) = xlx2 + xl + tanh(10(sin(5x2) - 2xi)). (4.6) 

The solution is anisotropic along the zigzag curve sin(5x2) — 2xi = and changes sharply 
in the direction normal to this curve(taken from |18tl22j). For more details see Table 4 
and Figure 5. 



Table 4: Four estimators and 5 in example 4.4 



step 


N 


E 




EI 


EIr 


\\H-Hr\\ 


5 


1 


146 


-0.644460 


0.106387 


0.842915 


0.082528 


520.885 




2 


325 


0.014871 


0.180764 


0.534598 


0.113211 


490.853 


0.15 


3 


756 


0.642053 


0.538530 


0.559720 


0.380420 


289.914 


1.25 


4 


1515 


0.940350 


0.854614 


0.935262 


0.766470 


105.952 


2.90 


5 


2826 


0.993081 


0.913258 


1.08005 


0.896111 


75.6472 


1.08 
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(a) 





Figure 5: (a) Four estimators, (b) the final mesh and (c) of the example 4.4. 
Example 4.5 This example is to solve the boundary value problem of 

- Au = 0, x(£n = (-0.5, 0.5) X (0, 0.5) U (-0.5, 0) x (-0.5, 0). (4.7) 

The Dirichlet boundary condition is chosen such that the exact solution is given by 

2 2 
u = r3 sm{ — 9), 
3 

where (r, 9) € Q are the usual polar coordinates. It is well known that the exact solution 
u G i73~^(r2)(Ve > 0). So we expect our estimators can be extended to more problems 
especially for those with low regularity. See Figure 6 for more details. 

Example 4.6 This example is to solve the boundary value problem 

-eAu + u = f, X € J] = (0, 1) X (0,1), (4.8) 

where the Dirichlet boundary condition and the right-hand side term are chosen such that 
the exact solution is the same as example 4.1 (taken from p9]). Note that 



e - UhWl^ + \\Viu - Uh)\\l^ 
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Figure 6: (a) Two estimators, (b) the final mesh and (c) of the example 4.5. 
(Er, EI and EIr are defined similarly). See Table 6 and Figure 7 for more details. 



Table 6: Four estimators and 5 in example 4.6 



step 


N 


E 


Er 


EI 


EIr 


\\H~Hr\\ 


5 


1 


94 


0.033119 


0.158037 


0.779738 


0.143494 


264.624 




2 


113 


0.049065 


0.216479 


0.811994 


0.206575 


229.060 


1.57 


3 


189 


0.666480 


0.568354 


0.584863 


0.448080 


148.367 


1.69 


4 


272 


0.958780 


0.867167 


0.918925 


0.756519 


58.9264 


5.07 


5 


278 


0.999348 


0.929424 


1.00642 


0.863337 


32.6764 


54.05 



From experiments above we conclude that our a posteriori error estimators r/ and r// are 
always asymptotically exact under various isotropic and anisotropic meshes. So we may 
guess that the superapproximation always holds during the adaptive process. 
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Figure 7: Four estimators of example 4.6. 



5 Conclusions 

In the previous sections we have developed a new type of a posteriori error estimators 
suitable for moving mesh methods under general meshes (especially anisotropic meshes). 
In our next paper we want to design adaptive algorithms using the estimators, i.e., to give 
a new metric tensor for moving mesh method. 

References 

[1] A. Agouzal, Y. Vassilevski, On a discrete Hessian recovery for PI finite elements, J. 
Numer. Math. 10 (2002) 1-12. 

[2] D. Ait-Ali-Yahia, W. Habashi, A. Tam, M.-G. Vallet, M. Fortin, A directionally 
adaptive methodology using an edge-based error estimate on quadrilateral grids, Int. 
J. Numer. Methods Fluids, 23 (1996) 673-690. 

[3] T. Apel, G. Lube, Anisotropic mesh refinement in stabilized Galerkin methods, Nu- 
mer. Math. 74(3) (1996) 261-282. 

[4] T. Apel, Anisotropic Finite Elements: Local Estimates and Applications, Advances 
in Numerical Mathematics, Stuttgart: Teubner, 1999. 

[5] R. E. Bank, R. K. Smith, Mesh smoothing using a posteriori error estimates, SIAM 
J. Numer. Anal., 34 (1997) 979-997. 

[6] R. E. Bank, J. Xu, Asymptotically exact a posteriori error estimators. I. Grids with 
superconvergence, SIAM J. Numer. Anal. 41(6) (2003) 2294-2312. 



14 



[7] R. Becker, An adaptive finite element method for the incompressible Navier-stokes 
equations on time-dependent domains, Ph.D. thesis, Ruprecht-Karls-Universitat Hei- 
delberg, 1995. 

[8] M. Berzins, A solution-based triangular and tetrahedral mesh quality indicator, SIAM 
J. Sci. Comput., 19 (1998) 2051-2060. 

[9] G. Buscaglia, E. Dari, Anisotropic mesh optimization and its application in adaptiv- 
ity, Internat. J. Numer. Methods Engrg. 40 (1997) 4119-4136. 

[10] G. Buscaglia, A. Agouzal, P. Ramirez, E. Dari, On Hessian recovery and anisotropic 
adaptivity. Fourth ECCOMAS Computational Fluid Dynamics Conference, Athens, 
1998, 403-407. 

[11] W. Cao, On the error of linear interpolation and the orientation, aspect ratio, and 
internal angles of a triangle, SIAM J. Numer. Anal. 43(1) (2005) 19-40. 

[12] M. J. Castro-Dfaz, F. Hecht, B. Mohammadi, O. Pironneau, Anisotropic unstructured 
mesh adaption for flow simulations, Internat. J. Numer. Methods Fluids 25(4) (1997) 
475-491. 

[13] L. Chen, P. Sun, J. Xu, Optimal anisotropic meshes for minimizing interpolation 
errors in L^-norm, Math. Comp. 76(257) (2007) 179-204. 

[14] P. G. Ciarlet, The finite element method for elliptic problems. Studies in Mathematics 
and its Applications, Vol. 4. North- Holland Publishing Co., Amsterdam-New York- 
Oxford, 1978. 

[15] L. Formaggia, S. Perotto, Anisotropic error estimates for elliptic problems, Numer. 
Math. 94(1) (2003) 67-92. 

[16] W. G. Habashi, M. Fortin, J. Dompierre, M. G. Vallet, Y. Bourgault, Anisotropic 
mesh adaptation: a step towards a mesh-independent and user- independent CFD, 
Barriers and challenges in computational fluid dynamics (Hampton, VA, 1996), 99- 
117, Kluwer Acad. Publ., Dordrecht, 1998. 

[17] F. Hecht, Bidimensional anisotropic mesh generator. Technical Report, INRIA, Roc- 
quencourt, 1997. 

[18] F. Hecht, A few snags in mesh adaptation loops. In proceedings of the 14th Interna- 
tional Meshing Roundtable, Springer- Verlag Berlin Heidelberg, 2005. 

[19] W. Huang, Metric tensors for anisotropic mesh generation, J. Comput. Phys. 204(2) 
(2005) 633-665. 

[20] G. Kuncrt, An a posteriori residual error estimator for the finite element method on 
anisotropic tetrahedral meshes, Numer. Math. 86(3) (2000) 471-490. 



15 



[21] Q. Lin, J. Lin, Finite Element Methods: Accuracy and Improvement, Science Press, 
Beijing, 2006. 

[22] K. Lipnikov, Y. Vasilcvski, Analysis of Hessian recovery methods for generating adap- 
tive meshes. In Proceedings of the 15th International Meshing Roundtable, Birming- 
ham, AL, September 2006, pages 163-171. 

[23] E. J. Nadler, Pieccwise linear approximation on triangulations of a planar region, 
Ph.D. Thesis, Division of Applied Mathematics, Brown University, Providence, RI, 
1985. 

[24] J. Pcraire, M. Vahdati, K. Morgan, O.C. Zienkiewicz, Adaptive remeshing for com- 
pressible flow computation, J. Comp. Phys. 72(2) (1987) 449-466. 

[25] M. Picasso, An anisotropic error indicator based on Zienkiewicz- Zhu error estimator: 
application to elliptic and parabolic problems, SIAM J. Sci. Comput 24(4) (2002) 
1328-1355. 

[26] J. R. Shewchuk, What Is a Good Linear Element? Interpolation, Conditioning, and 
Quality Measures, Eleventh International Meshing Roundtable (Ithaca, New York), 
pages 115-126, Sandia National Laboratories, September 2002. 

[27] K. G. Siebert, An a posteriori error estimator for anisotropic refinement, Numer. 
Math. 73(3) (1996) 373-398. 

[28] M. G. Vallet, C. M. Manole, J. Dompierre, S. Dufour, F. Guibault, Numerical com- 
parison of some Hessian recovery techniques. International Journal for Numerical 
Methods in Engineering, 72 (2007) 987-1007. 

[29] X. Zhang, Accuracy concern for Hessian metric. Internal Note, CERCA. 

[30] Z. Zhang, A. Naga, A new finite clement gradient recovery method: superconvergence 
property, SIAM J. Sci. Comput. 26(4) (2005) 1192-1213. 

[31] O. C. Zienkiewicz, J. Wu, Automatic directional refinement in adaptive analysis of 
compressible fiows, Intcrnat. J. Numer. Methods Engrg. 37 (1994) 2189-2210. 

[32] O. C. Zienkiewicz, J. Zhu, A simple error estimator and adaptive procedure for prac- 
tical engineering analysis. Int. J. Numer. Meth. Eng. 24 (1987) 337-357. 



16 



